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have included sufficient material to make this presentation self contained. Thus, all 
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other computer programs and used in different domains. Tables of numerical values 
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I. INTRODUCTION 

Some problems regarding the elliptical cylinders can be solved by using an analytical 
approach like that applied to circular cylinders: one separates the variables and the exact 
solution is given by expansions involving angular and radial Mathieu functions. These 
functions have been introduced by Emile Mathieu in 1868 by investigating the vibrating 
modes in an elliptic membrane [Ij. Details (tables or relations) concerning the Mathieu 
functions can be found for example in [2l|3llll|5l[6l[7l[8l|9l [TOl [TTl [12] . For circular cylinders 
the solutions involve readily available trigonometric and Bessel functions, while for elliptical 
cylinders there are still controversial and incomplete algorithms for computing the Mathieu 
functions. Largely applied computer programs provide only few or none routines refering to 
the Mathieu functions. 

A computational toolbox for Mathieu functions was implemented in Matlab p^. Since 
not all people are familiarized with the Matlab program, in this presentation the math- 
ematics is outlined. Tables of numerical values are provided. Following the outlines in 
this presentation, it would be a readily task to implement the Mathieu functions in other 
computer programs and use them in different domains. 

One reason for the lack of algorithms for Mathieu functions was probably the complicated 
and various notation existent in the literature. A main purpose for us was to simplify as much 
as possible the notation. With a simplified and self-contained notation, the use of Mathieu 
functions should be as simple as the use of Bessel functions. We largely followed the notations 
used by Stratton [6] and Stamnes [II],|12], but we introduced further simplifications. Since 
the notation and normalization for Mathieu functions vary in the literature, we have included 
sufficient material to make this presentation self contained. Thus, all formulas required to 
get the Mathieu functions are given explicitly. Tables of numerical values are provided. 
Examples of Mathieu functions applied to plane wave scattering by elliptical cylinders are 
given in [HI [15] . 
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II. FUNDAMENTALS 



A. Elliptical cylinder coordinates 



Let consider an ellipse in the plane {x,y) defined by equation (x/xq)^ + {y/yoY = 1 
with xq > yo. The semifocal distance / is given by p = — and the eccentricity is 
e = f/xQ < 1. The elliptic cylindrical coordinates {u,v,z) are defined by relations 



X = f cosh u cosv, y = f sinh u sin t>, 



z = z 



(1) 



with < M < oo and < f < 27r. In terms of {^,r],z), with ^ = cosh-u and t] = cosv, the 
elliptic cylindrical coordinates are defined by relations 



x = f^V, y = f^/{e-m-v'), 



z = z. 



(2) 



The contours of constant u are confocal ellipses (of semiaxes Xo = yo = /a/^^ — 1) and 
those of constant v are confocal hyperbolas. The z axis coincides with the cylinder axis. 
The scale factors hj, with j = C,,ti,z, are defined like as for any coordinate transformation 



h = f 



K = f 



1. 



(3) 



B. Wave equation in elliptic cylindrical coordinates 



The scalar wave equation (V^ + k^)U{r) = 0, where r is the position vector, k is the 
wave number, k = 27r-\/e/A, e is the permittivity, and A is the wavelength in vacuum, when 
expressed in elliptic cylindrical coordinates becomes 



/2 (cosh 2u — cos 2f) V dv"^/ dz"^ 



Using a solution of the form U = Z{z)S{v)R{u) gives 

/d^ 



U {u, V, z) 



0. 



(4) 



0, 



d^ 

df 2 

r d^ 



+ (a-2gcos2t;) S{v) = 0, 
- (a - 2g cosh 2m) 1 i?(M) = 0, 



(5) 
(6) 
(7) 
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where kz is the wave vector component on z direction, q = /c^/^/4, with = k'^ — k'^, and a 
is separation constant. Equation ^ has solution Z(z) = exp{ikzz). Equations ^ and ([T]) 
are known as the angular and radial Mathieu equations, respectively. 



III. ANGULAR MATHIEU FUNCTIONS 



In this presentation, only the periodic solutions of period vr or 27r are considered. For a 
given order n, there are four categories of periodic solutions satisfying (|6|, 

oo 

1 even-even: Seejv, q, n) = A^^J\q, n) cos{2jv), 

3=0 
oo 

2 even-odd: Seo{v, g, n) = J] AfJ+^\q, n) cos[(2j + l)v\, (8) 

oo 

3 odd-even: Soe{v^ q, n) = A'^^J'^ (g, n) sin(2jt>), 

i=i 

oo 

4 odd-odd: SUv, q,n) = Y, ^o'^'H?, n) sin[(2j + l)v]. 

j=0 

Apm with p,m = e,o are expansion coefficients. In the following, the angular Mathieu 
functions are denoted Spm{v,q,n), with p,m = e,o. Instead of two angular Mathieu func- 
tions, even Sep and odd Sop, with p = e,o pL2j, a single angular Mathieu function Spm, with 
p,m = e, o, is considered refering to all the four categories. For a given value of q there 
exist four infinite sequences of characteristic values (eigenvalues) a, for either value of a 
corresponding an infinite sequence (eigenvector) of expansion coefficients. 



A. Characteristic values and coefficients 

By subsituting ([s]) in (|6]), the following recurrence relations among the expansion coeffi- 
cients result 

1 even-even: 
aAi? - qA'^^J = 0, 



{a-A)A^^J -q[2A^^J +Aii^]=0 



\a 



(2jnA^^ - q[Afr'^ + Af^"-'^] = 0, J = 2, 3, 4 . . . (9) 



2 even-odd: 

(a - l)A^}] - q[A^^} + Af]] 



[a - (2j + lf]Afr'^ - q[A\T'> + A 



(2j-l) , .(2j+3)n 



0, J = 1,2,3- 



[a-(2jr]Ail^')-g[A^r'^ + A 



3 odd-even: 

(a - 4)Ai'e^ - gAit^ 

(2i~2) , .(2i+2)n 



J = 2,3,4. 



4 odd-odd: 

(a - + q[A^}} - Af}] 



[a - (2j + \f\A^^t^^ - q{A\r'^ + A 
recurrence relations can be written in matrix form [TT] . 



(2j-l) , .(2j+3)i 



1 even-even: 



( -a q 
2g 22 - a q 
q 42 - a g 00 
q e^-agO 

V ; ; ; ; ; ; 



J = l,2,3 



.(2) 
^ee 

4(4) 
^ee 

4(6) 



2 even-odd: 



/l 



+ g - a g 




















a q 

72 - a g 



/ A^ \ 



A3) 

.(5) 
.(7) 

V ^ I 



3 odd-even: 



( 1^ -a q 00 

g 42 - a g 00 

g & -a q 00 

g 82-agO 



/ Ail) \ 
4(4) 

.(6) 



oe 
(8) 



A 



6 





/l 

A. 


— n — n 


n 
H 


n 

W 


n 

w 


n 

w 






-^00 






Q 


6 — a 




n 
U 


n 
U 


n 

U ■ 




4(3) 

^00 


4 odd-odd: 







Q 


52 -a 


Q 





■ 




4(5) 












Q 


72- 


a q 


■ 




4(7) 

^00 




V 














••/ 


\ ■■ J 



0. 



(16) 



The matrices are real, tridiagonal, and symmetric for all categories, with the exception 
of the "1 even-even" category where the matrix is slightly non-symmetric. The eigenvalue 
problem is accurately solved in Matlab. In other computer programs it could be necessary to 
transform the slightly non-symmetric matrix in a symmetric one [TT|. Both the eigenvalues 
a and the corresponding eigenvectors {Apm, with p,m = e,o) are determined for either 
category at any order n. The order n takes different values for each category of Mathieu 
functions. For the purpose of avoiding any confusion, a distinction must be done between 
the n*^ order (in the succession of all orders) and the true value of that order. Thus, let 
denote n the order in the succession of all orders, and t the true value of order n. The values 
of n and t for the four categories of Mathieu functions are 



1 even-even: 

2 even-odd: 

3 odd-even: 
4 odd-odd: 



n = 0,1,2 



n = 0,1,2 



n = 1,2,3 



n = 0,1,2. 



t = 0,2,4- ■ 



t = 1,3,5- ■ 



it = 2,4,6- ■ 



t = 1,3,5-- 



Note that, if the notation is self-contained by all routines of Mathieu functions, there is no 
need to determine the specific values of n and t for either category of Mathieu functions 
since it is done automatically. 

B. Normalization and orthogonality 



Following [6l [TT] , the angular Mathieu functions are normalized by requiring that 

-dSop{v,q,n) 



Sep{0,q,n) = 1, 



dv 



v=0 



p = e,o. 



(17) 
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These requirements imply that, 

oo 

1 even-even: AfJ\q, n) — 1, 

3=0 
oo 

2 even-odd: J] ^) = 1> (18) 

i=o 

oo 

3 odd-even: ^ 2j7l(,2/)(g, n) = 1, 

oo 

4 odd-odd: ^(2j + n) = 1. 

i=o 

The orthogonahty relation for the angular Mathieu functions is 

/ Spm{v,q,n)Spm'{v,q,n)dv = Npm6mm', p,m,m' = e,o, (19) 

where Npm is normalization factor, 6mm' equals 1 if m = m' and equals otherwise. Then, 
the following relations for the normalization factor result, 

oo 

1 even-even: N^^iq, n) = 27r[4°)(g, n)]' + tt J2[Afj\q, n)]^ 

3=1 

oo 

2 even-odd: iVe„(g, n) = tt J^f^^^+^^lg, n)]^, (20) 

3=0 

oo 

3 odd-even: N^{q,n) = tt J][42^)(g, n)]^ 

oo 

4 odd-odd: Noo{q,n) - tt J][4^^+^)(g, n)]^ 

Since different normalization schemes have been adopted in the literature, much attention 
should be paid when numerical results provided by different authors are compared ones 
against the others. 



C. Correlation factors 

Let consider two regions of different permittivities, e and e'. The parameter q being 
different in the two regions, q ^ q', the characteristic values and expansion coefficients are 
also different. Let Spm and S' be the respective angular Mathieu functions. The correlation 
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factors Cpm{q,q',n), with p,m = e,o, between the angular Mathieu functions Spm and S'p^ 
are defined by relation 

/■27r 

Cp^{q,q',n) = Smm' Spm'iv,q,n)Sp^{v,q' ,n) dv, p,m,m' = e,o. (21) 

^0 

Using ([s]) gives 

1 even-even: C ee{q, q' ,n) = 2n A'^^J {q,n)A'J^^\q' ,n) 

oo 
oo 

2 even-odd: C^olg, g^ n) = vr ^ n), (22) 

CXD 

3 odd-even: CUq,q',n) = tt ^'oJKq , n) A'J~^^\q' , n) 

i=i 

oo 



4 odd-odd: C„,(g, g', n) = vr ^o'+'H?, n)A'i2^+i)(g', n). 

j=0 

D. Derivatives of angular Mathieu functions 

The derivatives of the angular Mathieu functions follow readily from ([s]) , 

1 even-even: — ^^^^ ' ^' — - = — 2jAfJ\q, n) sin(2jf ), 

2 even-odd: d^eo(^y,n) ^ _^(2j + l)42^+i)(g, n) sin[(2j + l)t;], (23) 



3 odd-even: '^^oe{v q.nj ^ J22jA^^j\q, n) cos(2j. 



oo 



dv 



4 odd-odd: ^^^ooiv^q^n) ^ j.^^^ ^ I) A'^^J+'\q,n) cos[i2j + l)v]. 

j=0 

IV. RADIAL MATHIEU FUNCTIONS 

Solutions of ([t]) can be obtained from ([s]) by replacing v by iu. Instead of sinf and cosw, 
the terms of the series now involve sinhw and coshw. The convergence is low unless 
is small. Better convergence of series results by expressing the solutions of ([T]) in terms of 
Bessel functions associated with the same expansion coefficients that are determined once for 
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both the angular and radial Mathieu functions. Either pair of angular and radial Mathieu 
functions are proportional to one another [B], 

Sep{iu, q, n) = V27rgep{q, n)Jep{u, q,n), p = e, o, (24) 

where Jgp are even radial Mathieu functions of the first kind and Qep are joining factors. 
When M = 0, 



5'ep(0,g,n) = 1, Jep{0,q,n) 
Thus, one obtains, 



1 



2ngep{q, n) 



p = e,o. 



(25) 



1 even-even: geeiq,n) = j^y- -5'ee(7r/2, g, ra), r = t/2, 

TTA),e{q,n) 



2 even-odd: geo{q,n) 



Similarly [6] 



■Ky/qAi^J{q,n} 



dSeo{v,q,n) 



dv J v=it/2 



r = {t- l)/2. (26) 



dJop{u, q, n) 



- iSop{iu,q,n) = V2TTgop{q,n)Jop{u,q,n), p = e,o. 

When u = 0, 

Jop{0,q,n) = 0, 
Thus, one obtains, 

3 odd-even: goe{q,n 



du iu=o y/27igop{q,n 



p = e,o. 



(27) 



(28) 



(2) 

7iqAle{q,n) 
(-1Y 

4 odd-odd: goo{q, n) = r- 5'oo(vr/2, q,n), r = {t - l)/2. 

'ny/qAlJ{q,n) 

Remember that t is the true value of order n. 



dv 



v=tt/2 



r = t/2. 



(29) 



A. Radial Mathieu functions of the first kind 



Since rapidly converging series are those expressed in terms of products of Bessel functions 
[TUt [TT] . in the following relations refer only to them. Similarly to the angular Mathieu 
functions, one may distinct four categories of radial Mathieu functions of the first kind 
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which are denoted Jpm{u, q, n), with p,m — e, o, 

1 even-even: Jgg{u,q,n) 

2 even-odd: Jeo{u,q,n) = a/| .fn/^' , ^(-1^^ ^Hg, 

+ J,(^;2)J,+iK)], r = (t-l)/2, (30) 

3 odd-even: J^eiu, q, n) = w| [~ ^ n)[J,_i(T;i) J,+i(t;2) 

- Jj-i(v2)Jj+i(vi)], r = t/2, 

4 odd-odd: Jo.Kg,n) = a/| .iV,^^^ , f;(-iy4^,^+^)(g, n)[J,(^i) J,-+i(^2) 



- J,(^2) r=(t-l)/2, 

where fi = ^cxp (— m) and ^2 = y^exp {u). The derivatives of the radial Mathieu func- 
tions of the first kind are 

1 even-even: r = t/2, 

-V2Jj{Vi)Jj+i{v2)], 

2 even-odd: r = (i-l)/2, 

- Jj+livi)Jj+i{v2)] + {2j + l)[Jj+i{vi)Jj{v2) - Jj{vi)Jj+i{v2)]Y 

3 odd-even: r = t/2, 
dJoe{u,q,n) [tt (-1) 



du 



a/I-L^ E(-iy+'4r'n?, ^)(4j + 4){ J,(^0^.-(^2) (31) 
cosh 2uJj+i (til) Jj-+i (112) - {j + l)[—Jj+i{vi)Jj{v2) + —Jj{vi)Jj+i{v2)]\, 
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4 odd-odd: r = (t - l)/2, 
dJooiu,q,n) fn (-1) 



+ J,+i(t;i)Jj+i(i;2)] - (2j + l)[J,+i(i;i)J,(?;2) + Jj{vi)Jj+i{v2)] \. 



B. Radial Mathieu functions of the second kind 

A second independent solution of ([t]) is obtained by replacing the Bessel functions of the 



first kind J„(f2) in (30) by the Bessel functions of the second kind Yn{v2) [IHlIilj- This 



solution is denoted Ypm{u, q, n), with p,m = e, o. 

1 even-even: YUu,q,n) = Y.{-iy Af^\q,n)J,{v,)Y,{v2). 

V Aee [q, nj 

r = t/2, 

2 even-odd: YUu,q,n) = .F- ^.{-1)' A^:j^'\q.n)[J,{v,)Y,^,{v2) 

+ r = (t-l)/2, (32) 

3 odd-even: YUu,q,n) = a/| t^^^ , f^i-lY Afj\q,n)[J,_,{v,)Y,^^{v2) 

~Y,^^{v2)Jj+i{vi)i r = t/2, 

4 odd-odd: YUu,q,n) = a/I t^^' , f^i'lY A^oo'^'\q,n)[J,{v,)Y,^,{v2) 

V ^A^^J{q,n) ^ 

-F,MJ,+i(t;i)], r = (t-l)/2, 
The derivatives of the radial Mathieu functions of the second kind are 
1 even-even: r = t/2, 

- V2JjiVi)Yj+i{v2)], 
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2 even-odd: r = (t-l)/2 
dYeo{u,q,n) Prr (-1) 



V ^ Alo'{q,n) >■ 



dw 

- J,+M)Yj+,{v,)] + (2j + l)[J,+iK)y,(t;2) - J,K)F,+i(t;2)] 
3 odd-even: r = t/2, 



' i=o 



cosh2uJj+i{vi)Yj+i{v2) - (j + Jj+i(t;i)Y^-(t;2) H Jj(i;i)Fj+i(t;2)] 

f 1 f2 



4 odd-odd: r = (t - l)/2, 
dFoo(M,g,?^ 

+ J,+MYj+,{v2)] - (2j + l)[J,^^iv,)Yjiv2) + J,K)F,+i(t;2)] k 



C. Radial Mathieu functions of the third and the fourth kinds 

Radial Mathieu functions of the third kind, analogous to the Hankel functions of the first 
kind are defined as follows P ITT] 

Hpmi{u,q,n) = Jpm{u,q,n) + iYpm{u,q,n), p,m = e,o. (34) 

Similarly, radial Mathieu functions of the fourth kind, analogous to the Hankel functions of 
the second kind are defined as follows [6l ITT] 

Hpra2iu,q,n) = Jpmiu,q,n) - iYpmiu,q,n), p,m = e,o. (35) 

V. IMPLEMENTATION OF MATHIEU FUNCTIONS IN MATLAB 

Following the notation of the four categories of angular Mathieu functions, the 
implementation in Matlab or in any other computer program is readily done by 
introducing a function code KF. The first step in any algorithm of Mathieu function 
computation is to find the characteristic values (eigenvalues) and the expansion coefficients 
(eigenvectors). In [13], this is done by routine "eig.Spm" which has q as input parameter 
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(see Table |T]). Besides g, the function code KF should be specified. Thus, if KF = 1, the 
routine "eig_Spm" solves the eigenvalue problem for category "1 even-even" of Mathieu 
functions, if KF = 2 for category "2 even-odd", and so on. The number of expansion 
coefficients is the same, it is set equal to 25, for all categories of Mathieu functions. 
Concerning the outputs of routine "eig_Spm" , fa is a line vector representing the 
characteristic values a for all the 25 orders; mc is 25 x 25 matrix, where the columns 
represent the eigenvectors (that is, the expansion coefficients) for all orders; vt is a column 
vector specifying the true value t for all orders. Note that the eigenvectors in mc were 



processed to obey equation (18). For the purpose to save the time of computation, all the 
other routines have mc as input (see Table |l]), the routine "eig_Spm" being called once, at 
the beginning of the computation, for any values of coordinates u and v that intervene in 
that computation. Since in many cases the convergence is assured by the first several 
orders, all the other routines have nmax < 25 as input. It means that those routines take 
into account only the first nmax orders, but for either order the length of the 
corresponding eigenvector is the same, equal to 25. The routine "extract _one .value" can be 
used to get a single value, and the routine "extract _one_column" to get a single 
eigenvector, corresponding to the order t. The derivatives of Spm, with p,m = e, o, are 
computed by routine "dSpm" . For both "Spm" and "dSpm" , v is expressed in radians, 
with values in interval (0,27r). The normalization, correlation, and joining factors are 
computed by routines "Npm" , "Cpm" , and "gpm" , respectively. The four kinds of radial 
Mathieu functions, Jpm, Ypm, Hp^i, and Hpm2, with p,m = e, o, are computed by routines 
"Jpm" , "Ypm" , "Hpml" , and "Hpm2" , respectively, and their derivatives with respect to u 
by routines "dJpm","dYpm","dHpml", and "dHpm2", respectively. 
Numerical values of the separation constant a, of the angular Mathieu functions Spm and 
their derivatives S'^^, with p,m = e, o, where the prime denotes differentiation with respect 
to f , are given in Tables |l![]lVl They can be compared with data in [2]. With the purpose 
to facilitate the comparison, since in [2] the normalization Npm = tt is applied, the data of 

are multiplied by ^/ir/Npm- 



Spm and S'p^ in Tables 



II 



IV 



Concerning the radial Mathieu functions, numerical values of Sep{iu,q,n) and 



-iSop{iu, q, n) are given for u = 0.5 in Tables V and VI They are multiplied by ^Jt^JN^ 



"pm 



and compared with data in P]. Note that Sep is correlated to the radial Mathieu function 



of the first kind Je„ by Eq. (24), whereas Sop is correlated to Jop by Eq. (27). We found 
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that, for parameters in [9], the values of Sep{iu,q,n) and —iSop{iu,q,n) calculated with 
Eqs. (24) and (27) differ from those obtained with Eq. dsl) by less than 7.5 x 10~^^. 
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TABLE I: Routines comprised in the toolbox |13| . 



Name of routine 



Routine call 



What the routine computes 



eig_Spm 



[va, mc, vt] =eig_Spm(i^F, q) 



Spm 

dSpm 

Npm 

Cpm 

Jpm 

dJpm 

gpm 

Ypm 

dYpm 
Hpml 

dHpml 
Hpm2 

dHpm2 

extract _one .column 
extract _one_value 



y=Spm{KF, v, mc, nmax) 
y=dSpm(KF, v, mc, nmax) 
y=Npm{K F, mc, nmax) 

y=Cpm{KF, mc, mc' , nmax) 

y=Jpm(KF, u, q, mc, nmax) 

y=A^pm{KF, u, q, mc, nmax) 
y=gpm{KF, q, mc, nmax) 

y=Ypm{K F, u, q, mc, nmax) 

y=dYpm{KF, u, q, mc, nmax) 
y=}ipml{KF, u, q, mc, nmax) 

y=dHpml(K F, u, q, mc, namax) 
y=Ilpm2(K F, u, q, mc, nmax) 

y=dRpm2{K F, u, q, mc, namax) 
y=extract_one_column(i^i^, t, mc) 
y=extract_one_value(isrF, t, vec) 



Vector of characteristic values va, matrix of 
coefficients mc, and vector of orders vt, at given 
function code KF and elliptical parameter q >0. 
Angular Mathieu functions Spm, [Eq. (8)]. 



Derivatives with respect to v of Spm, [Eq. (23)] 



Normalizing factors of angular Mathieu functions 



Spm, [Eqs. (19) and (20)]. 



Correlation factors of Spm and 5"^^, having matrices 



of coefficients mc and mc' , [Eqs. (21) and (22)]. 



Radial Mathieu functions of the first kind J, 



pm. 



[Eq. (30)] 



Derivatives with respect to u of Jpm, [Eq. (31)] 



Joining factors for pairs of angular, Spm and 



radial, Jpm Mathieu functions, [Eqs. (24)-(29)] 



Radial Mathieu functions of the second kind 



Ypm, [Eq. (32)]. 



Derivatives with respect to u of Ypm, [Eq. (33)]. 



Radial Mathieu functions of the third kind Hpmi , 



[Eq. (34)] 



Derivatives with respect to u of Hpmi ■ 

Radial Mathieu functions of the fourth kind Hpm2, 



[Eq. (|35[]. 

Derivatives with respect to u of Hpm2- 
Extracts one column from mc at given t. 
Extracts one value from vec at given t. 



TABLE II: Values of See multiplied by 7ee = \/ tt /Nee to be compared with data 



J. 

t 


q 


a 


7ee5'ee(0,g',n) 


7eeOee(7r/2,q',nj 


u 


u 





0.7071067811865 


0.7071067811865 




5 


f o r\r\r\ A r\(^r\o r" '\ f 

-o.8000460208olo 


0.0448001816519 


1 ^^AOAor'''7Ar'f\or\ 

1.3348486746980 




1 n 
iU 


-io.yoDy /yyoooooy 


u.uu/DZDol 


1 /I '7n'7i on 
i.4DoODU4 / U / izy 




io 


-zz.oloUoi /DUoD4U 


U.UUlyozoUooloz 


1.00U1Uo14DDoDD 




on 


-oi.olooyUU ^ UooD4 


U.UUUbUo ( 4oozyz 


i.DUyoyUoo { oyoy 




oc; 
zo 


/in O l^fi'7'7n 1^ /I 
-4U.Z0D/ /yo400DO/ 


U.UUUzi0o0oUlo4 


l.DO /OlUzyooZOO 


1 


u 


A r\c\c\c\c\c\r\c\r\c\c\r\c\ 

4.0000000000000 


1.0000000000000 


1 nnnnnnnnnnnnn 

-1 .0000000000000 




5 


7.4491097395292 


n '70 r on /I ono /I nnr^ 

0.7352943084007 


-0.7244881519677 




1 n 
iU 


(.ill od9o49/ /9d 


0.z4oooo349z913 


-0.9zD /59zD4izD3 




1 ^ 

io 


O.U/ /^9ooi9 /^04oo 


U.U to /^9zoz ( o4d4 


-i.Ui99DDZZDUoUo 




on 


i.io4zozooi3z4Do 


U.UzoD4o943i4 ( i 


-i.U/0z9ozzo/ 




zo 


-O.0ZZiD4 i Zl lOOO 


U.UiiuizoDDooUy 


-i.iiDz /oy Oozy Oo 




U 


100.0000000000000 


1 nnnnnnnnnnnnn 
1.0000000000000 


1 nnnnnnnnnnnnn 
-1 .0000000000000 







iUU.izDoD9ziDiDob 


1 nor^nn^novnoo/i 
i.Uzo99oUz /Uo94 


-U.y / 034 / 4o / z3bU 




10 


100.5067700246816 


1.0538159921009 


-0.9516453181790 




15 


101.1452034473016 


1.0841063118392 


-0.9285480638845 




20 


102.0489160244372 


1.1177886312594 


-0.9057107845941 




25 


103.2302048044949 


1.1562399186322 


-0.8826919105637 



TABLE III: Values of Seo and 3'^,^ multiplied by 7eo = \/ vr /Neo to be compared with data 



t 


q 


a 


7eoSeo{0,q,n) 


7eo<3eo(,7r/2,g,nj 


1 


u 


1 nnnnnnnnnnnnn 

1.0000000000000 


1 nnnnnnnnnnnnn 

1 .0000000000000 


1 nnnnnnnnnnnnn 

-1 .0000000000000 




5 


1.8581875415478 


0.2565428793224 


o /i/^nn^ionno/inr'^ 

-3.4690420034057 




1 n 
iU 




U.Uoooyo/4/ /4/Z 


A Qt^r\AQQ'iC\AACiCiA 

-4.ooU4oooU44yo4 






-o.iUiiU0ioiD4io 


U.U10U4UU0D40oo 


-0. ( 04zUD4oyU0lU 




on 


-14.4yioUi4zoi ( 4o 


U.UUoUoioio /d4 ( 


-D.4yU000 (oZoovv 




oc; 
zo 


-zi.oi4oyyoyuDDo ( 


u.uuiyiiuoioUD ( 


- / .iUb /4izoozyUi 


5 


(J 


o c! r\r\r\r\r\r\r\r\r\r\r\r\r\ 

25.0000000000000 


1 nnnnnnnnnnnnn 

1 .0000000000000 


r; nnnnnnnnnnnnn 

-5.0000000000000 




5 


25.5499717499816 


1.1248072506385 


-5.3924861549882 




1 n 
iU 


z / . (Uo /bo 1 oo9o9o 


i.z0oUi994ioUoo 


Q01 ovftc^/11 1 ftnn 
-O.oziz /bu4iibUy 




io 


Qi nc^vQoi oc^oi von 
oi.9i3 ioZLZoZi. 


1 1 OQ/l QOOQn/l 1 Q1 

i.i9o4ozzoU4ioi 


1 1 ni /I no QQ /I nft /I 
-0.iiyi4yooo4Ub4 




on 


Qfi ft /I /I nQO 1 QOQ 

od.d449o9 i o4ioz0 


U.90D0 ( OOoi4zzb 


c; WQfivcoc^nnft/1 /I 
-0./ /^ob/^t)Z0UUb44 




zo 


/in nc^m nnnc ^Qnw 


n fii nfiO/iQi nnr^nv 
U.DiUDy4oiUUuU ( 


V nr^OQQ/lOOl ftc;c:Q 


lb 


U 


ooc^ nnnnnnnnnnnnn 
225.0000000000000 


1 nnnnnnnnnnnnn 
1 .0000000000000 


1 nnnnnnnnnnnnn 
15.0000000000000 




r 



oor^ n^^Qi O/l Vfivnoft 
zzo.UOOoiz4 /b/U9D 


1 ni 1 ooQ'zQor^onfi 


1 1 fiQfir^V/l vonftno 
io.ibc)bo (4/ zUbUz 




10 


225.2233569749644 


1.0228782824382 


15.3198803056623 




15 


225.5029562446541 


1.0347936522369 


15.4687435032830 




20 


225.8951534162079 


1.0470843441629 


15.6102785232380 




25 


226.4007200447481 


1.0598004418139 


15.7444725050679 



TABLE IV: Values of Sop and S'^^ multiplied by ^ir/Nop^p = e,o (see p]) 













t 


q 


a 


y/TT/NoeS„^{0,q,n) 


/ /AT 0/ / /n \ 


2 





4.0000000000000 


2.0000000000000 


e~\ in in in in /n /n /n /n /n /n /n in in 

-2.0000000000000 




5 


2.0994604454867 


0.7331661960372 


-3.6405178524082 




10 


o ooo"i r" oo o rri f '7r\ 

-2.3821582359570 


0.2488228403985 


-4.8634220691653 




15 


-8.0993467988959 


0.0918197143696 


-5.7655737717278 




20 


-14.4910632559807 


0.0370277776852 


^n /n r' (n (n (n ^i/no'^o 

-6.4907522240373 




25 


-21.3148606222498 


0.0160562170491 


-7.1067719073739 


10 





100.0000000000000 


10.0000000000000 


-10.0000000000000 




5 


100.1263692156019 


9.7341731518695 


-10.2396462566908 




10 


100.5067694628784 


9.4404054347686 


-iin Ar't~tt~\A'~7r't~t't^ACir'' 

-10.4539475316485 




1 K 
10 


lUi.i4oi i ZZ\)Z\){J\)Z 




in ft AOQCiC\Q'7'7fit^^Q 

-lU.D4z8yyo/ (bobo 




on 




0. /ooo4oUoUlooU 


in on '70 /1 1 '70 1 ooc^ 




25 


103.2256800423735 




-i(j.y4i<35oMUoiyi 












t 


q 


a 


^/tt/NooS^^{0, q, n) 


/ /AT / /O \ 

^/TT/NooSoo[TT/2,q,n) 


1 





1.0000000000000 


1.0000000000000 


1 in /n in /n /n in /n in in /n in in /n 

1.0000000000000 




5 


-5.7900805986378 


0.1746754006198 


1.3374338870223 




10 


-13.9365524792501 


0.0440225659111 


1.4687556641029 




15 


-22.5130034974235 


0.0139251347875 


1 r'rT\'i 1 r" A r" ^ f* 

1.5501150743576 




20 


-31.3133861669129 


in in r' in ii^ninin-i 

0.0050778849001 


1.6098915926038 




25 


-40.2567789846842 


0.0020443593656 


1.6575103983745 


5 





25.0000000000000 


r" /n in in /n in in /n in in in in in in 

5.0000000000000 


1 in in in in in in in in in in in in in 

1.0000000000000 




5 


25.5108160463032 


oo/nr'^T/n/n-1 /n ^*n ^/^ 

4.3395700104946 


/n *n /n /n ^n in (n in 4n /I 

0.9060779302024 




10 


26.7664263604801 


3.4072267604013 


0.8460384335355 




15 


27.9678805967175 


2.4116664728002 


/n oo'^^n^/nOiiin/n'iinr' 

0.8379493400125 




20 


28.4682213251027 


1.5688968684857 


0.8635431218534 




25 


28.0627658994543 


0.9640716219024 


0.8992683245108 


15 





225.0000000000000 


1 r' in in in in in in in in in in in in in 

15.0000000000000 


1 in /n /n in /n in in /n in /n in in in 

-1.0000000000000 




5 


225.0558124767096 


14.8287889732852 


-0.9889607027406 




10 


225.2233569749643 


14.6498600449581 


-0.9781423471832 




15 


225.5029562446537 


14.4630006940372 


-0.9675137031855 




20 


225.8951534161767 


14.2679460909928 


-0.9570452540613 




25 


226.4007200438825 


14.0643732956172 


-0.9467086958781 
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TABLE V: Values of Sep{iu, q, n) for u = 0.5 multiplied by y^vr/A^'ep, where p = e,o, compared 
with data in [9^ 





q 


Values at p — 6 


DaAjo. in [y 


-f 


q 


Values at p — o 


uaTa lu. [yj 


n 
u 





n ni Q'io^'i^\AQ^ nn7i 


n ni QQ9 
-u.uiyoz 


i 





U.UZ144U ( 4oloOOZ ( 


n noi A A 

U.UZi44 




iU 


-U.UU lUooZov I ioiyo 


-U.UU / Uo 




1 n 


n riQQfiQ/l OQ7/i P^QP^OP^ 

-U.UoODo4zo 1 4oo0zo 


-U.UooDo 




on 


-U.UUUiDy4ii4iO 1 oO 


n nnni 
-U.UUUlO 




on 


-U.UUoo 1 oooooUyD4z 


n nnQQ7 
-U.UUoo ( 


o 
z 


t3 


U.44Dyo 1 400 / 4iUDo 


Ci A A fiOQ 
U.44oyo 


Q 

o 





i.zUOOzozD / UODooo 


i.zUoo 




10 


-0.063855921612085 


-0.06385 




10 


0.235940782144547 


0.23594 




20 


-0.024916657795101 


-0.02491 




20 


-0.097385461808731 


-0.09738 


4 


5 


2.234088244534832 


2.2341 


5 


5 


3.864089377116713 


3.8641 




10 


1.039103163573830 


1.0391 




10 


2.285610444240526 


2.2856 




20 


-0.143991090269732 


-0.14399 




20 


0.274270780278172 


0.27427 



TABLE VL Values of —iSop{iu, q, n) for u = 0.5 multiplied by ^JiiJWfyp^ with p = e^o, compared 
with data in |9] 



t 


q 


Values at p = e 


Data in [9^ 


t 


q 


Values at p = o 


Data in [9 


2 


5 


0.238342768735937 


0.23834 


1 


5 


0.036613617783886 


0.03661 




10 


0.028675814044625 


0.02867 




10 


0.000750806874015 


0.00075 




20 


-0.003176296415956 


-0.00317 




20 


-0.000538258353937 


-0.00053 


4 


5 


1.883560277440876 


1.8836 


3 


5 


0.806555153528872 


0.80655 




10 


0.769679129538722 


0.76968 




10 


0.204495885546638 


0.20449 




20 


0.040515136278697 


0.04051 




20 


-0.005279473480675 


-0.00527 


6 


5 


6.6066602369876 


6.6067 


5 


5 


3.667530204538722 


3.6675 




10 


4.1161420952367 


4.1161 




10 


1.972361938552091 


1.9724 




20 


1.1805904286267 


1.1806 




20 


0.320398855944192 


0.32040 



